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^- 1 Abstract 

Q> ■ The effect of the anisotropy in the exchange interaction on the phase diagrams 

and magnetization behavior of the Heisenberg thin film has been investigated 
with effective field formulation in a two spin cluster using the decoupling ap- 
r^ \ proximation. Phase diagrams and magnetization behaviors have been obtained 

O ■ for several different cases, by grouping the systems in accordance with, whether 

a^ \ the surfaces/interior of the film have anisotropic exchange interaction or not. 

Particular attention has been paid on the evolution of the special point coor- 
dinate with the anisotropy in the exchange interaction. It has been concluded 
that for increasing ratio of the anisotropics in the exchange interaction between 
t/3 I the surface and interior of the film, special point disappears at a special value 

of this ratio. 
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rj^ ■ 2 Introduction 

Q . Recently, there has been growing interest both theoretically and experimentally 

O ' on the finite magnetic materials especially on semi-infinite systems and thin 

films. The magnetic properties of the materials in the presence of the free 
surfaces are drastically different from the bulk counterparts. This is because 
^ . of the fact that, free surface breaks the translational symmetry, i.e. surface 

Cn ' atoms are embedded in an environment of lower symmetry than that of the 

inner atoms [3 [5] . If the strength of the surface exchange interaction is greater 
than a critical value, the surface region can exhibit an ordered phase even if the 
bulk is paramagnetic and it has a transition temperature higher than the bulk 
one. This fact has been observed experimentally[3l|4j[5]. 
I^ . The development of the molecular beam epitaxy technique and its applica- 

tion to the growth of thin metallic films has stimulated renewed interest in thin 
film magnetism. In a thin film geometry, one important question is dependence 
of the thermodynamic properties and Curie temperature of the film on the thick- 
^^ , ness of the film. It was experimentally found that, the Curie temperature and 

H ' the average magnetic moment per atom increases with the increasing thickness 

- - - of the film [51 [7] . This thickness dependent Curie temperature also has been 

measured in Co[H], Fe[n] and Ni[Tn] films. 

In order to get insight on the phase transition characteristics and thermo- 
dynamic properties of thin films, it is valuable in theoretical manner to solve 
related models in thin film geometry such as Hesinberg and Ising models. One 
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class of the films which exhibit a strong uniaxial anisotropy[TT] can be modeled 
by Ising model. These systems have been widely studied in literature by means 
of several theoretical methods such as Monte Carlo (MC) simulations [H] , mean 
field approximation (MFA) [T3j and effective field theory (EFT) [Hj. In the Ising 
case, the study of ferroelectric films treated by transverse Ising model (TIM) 
are more common in the literature since fabricated ferroelectric thin films such 
as Barium titanatc BaTiOz and Strontium titanate SrTiO^ [15l [161 HI] have 
a lot of technological applications in optoelectronics [IH] and microelectronics 
[T5] . The TIM on the thin film geometry can be solved by a variety of tech- 
niques such as MFA [ID] and EFT HI]. Besides, within the EFT formulation, 
the effect of the dilution on the phase diagrams and magnetic properties of the 
ferroelectric thin films (by means of the TIM) has been studied [22] and also long 
range interactions have been taken into account |23| . In order to simulate more 
realistic models, different exchange interactions have been defined and different 
transverse fields have been considered on the surface and bulk, which anticipates 
the mimics of the surface effects within the framework of EFT [53]. Another 
improvement of the Ising thin film models is the consideration of amorphisation 
of the surface due to environmental effects, and these realizations have also been 
solved using EFT [5S]. There are also higher spin Ising thin films, e.g. spin-1 
Ising thin films have been studied [13 HZ] ■ 

Although the Ising model has been studied on the thin film geometry widely, 
less attention was paid on the Heisenberg model on thin film geometry, in com- 
parison with the Ising model. Thin films which do not exhibit a strong uniaxial 
anisotropy requires to solve the Heisenberg model in the thin film geometry. In 
order to see the effect of the presence of the surface in the system on the critical 
and thermodynamical properties, Heisenberg model on a semi infinite geome- 
try has been solved using a wide variety of techniques such as Green function 
method [28], renormahzation group technique [29], MFA [S^, EFT [3n[32ll33] . 
high temperature series expansion [34] |35]. As in the case of Ising counterpart, 
Heisenberg model in a thin film geometry has been solved in a limited case. For 
instance, Heisenberg model on a thin film geometry with Green function method 
[MJISIIISH], renormahzation group technique [33, MFA [30], EFT [HSI] and 
MC [Hm, are among them. 

The aim of this work is the determine the effect of the anisotropy in the 
exchange interaction on the phase diagrams and magnetization behavior of the 
Heisenberg thin film. For this aim, the paper is organized as follows: In Sec. [3] 
we briefly present the model and formulation. The results and discussions are 
presented in Sec. U] and finally Sec. |S] contains our conclusions. 

3 Model and Formulation 

Thin film can be treated as a layered structure which consist of interacting 
L parallel layers. Each layer is defined as a regular lattice with coordination 
number z. When we choose z = 4, this means that each layer has square lattices 
and each nearest neighbor layer have interaction. The Hamiltonian of the thin 



film is given by 

^ = - E (4<^^ + 4^^" + 4^?^) - E^^? (1) 

where sf,s^ and sf denote the Pauh spin operators at a site i. Jfj^Jfj and 
Jf- stand for the anisotropy in the exchange interactions between the nearest 
neighbor spins located at sites i and j and Hi is the longitudinal magnetic field 
at a site i. The first sum is carried over the nearest neighbors of the lattice, 
while the second one is over all the lattice sites. The exchange interaction 
components {Jfj^ Jf^, Jf,) between the spins on the sites i and j take the values 
according to the positions of the nearest neighbor spins. The two surfaces of 
the film have the intralayer coupling components (Jf, Jj*, Jf). The interlayer 
coupling between the surface and its adjacent layer (i.e. layers 1, 2 and L — 1, L) 
is denoted by {J2 ,J2,J2)- For the rest of the layers, the interlayer and the 
intralayer couplings are assumed as {J2,J2,J2)- 

We use the two spin cluster approximation as an EFT formulation, namely 
EFT-2 formulation |44| . In this approximation, we choose two spins (namely si 
and S2) in each layer and treat interactions exactly in this two spin cluster. In 
order to avoid some mathematical difficulties, we replace the perimeter spins of 
the two spin cluster by Ising spins (axial approximation) |45j . After all, by using 
the differential operator technique and decoupling approximation (DA) |46| . we 
can get an expression for the magnetization per spin, i.e. Jn = (^ {sf + Sj)). In 
the thin film geometry number of L different representative magnetizations for 
the system by following the procedure given in Ref. |41| can be given as, 

TOi = {Qi^i02.2) .flix,y,Hi,H2)\x=0,y=0 

rrik = (Q2,k-i'c>l^k'^2.k+i) f2{x,y,Hi,H2)\x=o,y=o,k^2,3,...,L-l 

■mi = (02,L-10?l)/i (2^,2/,ii^l,i^2)U=0,y=0■ 
(2) 

Here rrii, {i — 1,2, ... ,L) denotes the magnetization of the i*'* layer. The oper- 
ators in Eq. ([2|) are defined via 

Ok,i = [Akx + miBkx] [Aky + miBky] (3) 

where 



Akm = cosh(J|Vm) 

Bkm = sinh(J|Vm), k ^ l,2;m = x,y. 

The functions in Eq. ^ are given by 



(4) 



fn{.x,y,Hi,H2) = 



x + y + Hi+H2 sinh (/?X^"^ 

^,$"' cosh (/?X^"^) + exp (-2^ J-) cosh (/3yJ"^) 

(5) 



where 



^0 



{J:^-jyr + {x + y + H,+H2f '^' 



iJ^ + jyf + {x-y + Hi-H2 



1/2 



(6) 



with the values n = 1,2. In Eq. ([S]), we set /3 = 1/(A:bT) where ks is Boltzmann 
constant and T is the temperature. 

Magnetization expressions given in closed form in Eq. ^ can be constructed 
via acting differential operators on related functions. The effect of the exponen- 
tial differential operator to an arbitrary function F{x) is given by 



exp (aV)F (x) ^ F {x + a) 



(7) 



with any constant a. 

With the help of the Binomial expansion, Eq. ^ can be written in the form 
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*i = EEEE^i(P'9'^'^)"^r'^2+' 



p=0g=0r=0s=0 
113 3 11 
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EEEEEE^2 (P, q, r, s, t, v) mt\ml+^mlX\ (8) 



p=0g=0r=0s=0t=0'U=0 
3 3 11 



•m-L = EEEE-^''^! (^' 5' ''' *) ^^L 
p=0g=0r=0s=0 
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where 



Ki {p, q, r, s) 
K2{p,q,r,s,t,v) 




3-p A3-q ,1-r Al-snP 



At7K M^''A2y'Bf^BlyB2.Blyh ix,y,H,,H,) U=o,y=o 



5-(p+r+t) .4-{q+s+v) p+r+t r^q+s+v 
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'f2ix,y,Hi,H2)\x=0,y=0- 

(9) 



These coefficients can be calculated from definitions given in Eq. Q by using 
Eqs. © and (HI). 

For a given Hamiltonian parameters and temperature, by determining the 
coefficients from Eq. (j9]) we can obtain a system of coupled non linear equations 
from Eq. ([8]), and by solving this system we can get the longitudinal magneti- 
zations of each layer {rrii, i = 1, 2, . . . , i). The total longitudinal magnetization 
(m) can be calculated via 



1=1 



(10) 



Since all longitudinal magnetizations are close to zero in the vicinity of the 
second order critical point, we can obtain another coupled equation system to 
determine the transition temperature by linearizing the equation system given 
in Eq. ^, i.e. 



771 1 



nik 



rriL 



= [1^1(1,0, 0,0) + 7^1 (0,1, 0,0)] 7711 + 

[ifl(0,0,l,0)+i^l (0,0, 0,1)] 7772 

= [if3(l,0,0,0,0,0) + ii'3(0,l,0,0,0,0)]777fc_i- 
[Ks (0, 0, 1, 0, 0, 0) + K3 (0, 0, 0, 1, 0, 0)] mk+ 
[K3 (0, 0, 0, 0, 1, 0) + K3 (0, 0, 0, 0, 0, 1)] nik+i 

= [ifi(l,0,0,0) + /vi(0,l,0,0)]77iL + 

[Kl(0,0,l,0)+/Vl(0,0,0,l)]777L_l. 



(11) 



Critical temperature, as well as critical transverse field can be determined 
from det(A) = where A is the matrix of coefhcients of the linear equation 
system given in Eq. pT|) . 

4 Results and Discussion 

In this section, we treat the anisotropic Heisenberg model on the thin film 
geometry for several cases, with no magnetic field (i.e. iJi = iJ2 = in the 
formulation). The system has two different exchange interactions which control 
the spin-spin interactions between the spins belonging to interior (^fi'^l'i'^l) 
and the surfaces ( Jf , Jf^Jf) of the film. According to the relations between the 
components of these exchange interactions, let us entitle the case Jf = Jf = Jf 
and (Jf = J| = J2) as a thin film with isotropic surfaces (interior), and the 
case Jf = Jf, Jf = 1 and (Jf = J2,Ji = 1) as a thin film with XXZ type 
surfaces (interior). In other words, we group the portions of the film (interior 
and surfaces) according to the symmetry properties of the relevant exchange 
interactions as isotropic or XXZ type. 

Let us select the unit of energy as J ( J > 0) and scale all components of the 
exchange interaction with J, 

< = y, i = l,2,z/ = x,|/,z (12) 

here since we deal with the ferromagnetic interactions, defined dimcnsionless pa- 
rameters in Eq. (|12|) arc all positive or zero. Let the parameter g^, (1/ = x, y, z) 
controls the ratio of the same components of exchange interactions between the 
surfaces and other layers, 

<li' = -^, '^^x,y,z (13) 

'"2 

where r^ stands for the scaled (dimcnsionless) v component of the exchange 
interaction between the spins, which belong to surface, while r2 represents the 
V component of the exchange interaction between the remaining spins. 

4.1 Isotropic Model 

In this case, all components of the exchange interactions are equal to each other, 
i.e. 

n. 1 = 1,2 (14) 



and choose the unit of energy such that r2 ~ 1, i.e. unit of energy in Eq. (|12[) 
is J = J|(= Jf = Jf) . The film consist of isotropic surfaces on an isotropic 
interior. 

The thermodynamic properties of the system are affected by three param- 
eters, which arethe film thickness (L), ratio between the surface intralayer ex- 
change interaction and the other exchange interaction (ri), and the temperature 



{ksT/ J). Thus the critical temperature of the system depends on the param- 
eters L and ri. In a system with a surface, a typical phenomena occurs and 
this phenomena shows itself in a thin film geometry. The phase diagrams for 
different film thickness (L) intersect at a special point which can be denoted by 
(/cbT*/J, r*) in the (fcsTc/J, ri) plane. For the values of ri > rl, thin films 
have higher critical temperature than the corresponding bulk system while the 
reverse is valid for ri < r\. This means that magnetically disordered surface 
can coexist with a magnetically ordered bulk phase for the values of ri that 
provide ri < r\ while for the values ri > r\, surface can reach the magnetically 
ordered phase before the bulk. 

For the value of ri = r\, thin film with thickness L has the same critical 
values as that of the bulk, independent of the L. The value of ri also changes 
the relation between the film thickness and the critical temperature of the film. 
The thicker films have higher critical values for ri < rj while thicker films have 
lower critical values for ri > r* than the thinner ones. 

The phase diagrams for the isotropic Hciscnbcrg model in a thin film geome- 
try in a {ksTc/J, ri) plane can be seen in Fig. ([T]). Above mentioned typical be- 
havior can be seen in this figure. For the special point coordinate {ksT* /J,rl), 
EFT-2 formulation gives (4.8910, 1.3454), where the first one is just the critical 
temperature of the corresponding bulk system (the system with simple cubic 
lattice) in the same model[ll]. This value can be compared with the same 
system within the high temperature series expansion method (4.395, 1.325) [3S] 
and the value of the classical counterpart of the same system within the same 
formulation (5.030, 1.332) [4^. It is a known fact that the classical system has 
slightly higher critical temperature. Also, comparison of the values shows that 
EFT-2 formulation gives higher value for the ksT* / J and r*, according to the 
high temperature series expansion method. 

One of the efficient way of describing the behavior of the order parameter 
with temperature and ri is to treat the equally valued magnetization curves 
in the (fc^T'c/J, ri) plane. In Fig. ©, we depict this for the film thickness 
values of L = 3 and L = 10. As seen in Figs. ^ (a) and (b) the variation 
of the magnetization with these two parameters is qualitatively the same. But 
the curves become almost parallel lines with the ri axes when film thickness 
increases. This means that, magnetization value gets more and more affected 
by ri, when the film thickness decreases. 

The relation between the critical temperature of the film and the correspond- 
ing bulk system shows itself also in the relation between the magnetization of 
the surface layer and middle layer. In other words, chosen ri also determines 
the relation between the magnetization of the surface layer and middle layer 
for any temperature, which is below the critical temperature. As an example 
of this situation, the variation of the magnetization of the surface and middle 
layer as well as the film with temperature can be seen in Fig. ([3]) for L = 7, 
(a) ri = 0.5 < rl and (b) ri = 2.0 > r*. As seen in Fig. (O (a), for the values 
of ri < rj , surface layer has lower magnetization than the middle layer and the 
situation gets reverse when ri > r^ (Fig. ([3]) (b)). Of course this situation is 
valid for the temperatures that are not so close to zero. Even though ri = 0.0, 



the ground state of the surface layer will be ordered and the value is 1.0. In 
the absence of the interaction between the spins on the surface layer (ri = 0.0), 
the surface is completely ordered due to the interaction of the surface with the 
neighbor layer. The order of the inner layers dictates the order of the surface 
layer in the isotropic interaction case. 
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Figure 1: Variation of the critical temperature with r\ for some selected film 
thickness values of L = 3,4,6,10 and r^ = 1.0. The horizontal line represents 
the critical temperature of the isotropic Heisenberg model with the simple cubic 
lattice structure. 




Figure 2: Equally valued magnetization curves of the isotropic Heisenberg model 
in a (fcsT/J — ri) plane with the film thickness values (a) L = 3 and (b) L = 10 
and r2 = 1.0. The values of the magnetizations have been underlined near the 
related curve. 



4.2 Anisotropic Model 



In this case the anisotropics in the exchange interactions are related to each 
other as 

r-=rf = n, r^ = l, * = 1,2 (15) 
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Figure 3: Variation of the magnetization of the surface layer (nii), inner layer 
(7714) and total (m) magnetization with temperature for L = 7 and selected 
values of (a) ri — 0.5 and (b) j'l = 2.0 and r2 = 1.0. 



In order to investigate the effect of the anisotropy in the exchange interaction 
on the critical properties of the system, we choose for the z components of the 
exchange interactions as rf = r| = 1. This means that we choose the unit of 
energy in Eq. (fT2|) as J = J|(= Jf). From Eq. ([T3| and according to Eq. (fTSJ) . 
Qx = Qy and let q — Qx — Qy Thus, the parameter q controls the ratio of the 
anisotropies between the surfaces and other layers. We can treat this case in 
two sub-cases whether the interior of the film is isotropic or not. 



4.2.1 XXZ-type surfaces on an isotropic interior 

In this case, the inner layers have isotropic exchange interaction both along the 
interlayer and intralayer and, in order to focus on the effect of the anisotropy in 
the exchange interaction which comes from the surfaces, let us choose ^'2 = ^2 = 
1. In other words, the film consist of a body which has isotropically interacting 
spins, coated with surfaces containing intralayer spin-spin interactions which are 
XXZ-type symmetric. Thus, the anisotropy in the exchange interaction of the 
system comes only from the surfaces of the film. In order to see the effect of the 
anisotropy in the exchange interaction which comes from the surfaces, we depict 
the phase diagrams in a (fc^Tc/ J, q) plane in Fig. (|4]). At first sight, we can see 
from Fig. (J4]) that, rising anisotropy in the exchange interaction decreases the 
critical temperature of the film. Most affected films are the thinner ones (com- 
pare the curves labeled by A and D in Fig. (|4])). Since the surface and the inner 
layers are interacting via the exchange interaction J2, the effect of anisotropy 
can diffuse the inner layers of the film and effect the critical temperature of the 
system. This propagation becomes easier for the thinner films in comparison 
with the thicker ones, as seen in Fig. (j4]). Also, as shown in Fig. (j4|), all crit- 
ical temperatures lie below the critical temperature of the corresponding bulk 
system of simple cubic geometry and it is shown as a thinner horizontal line in 
the diagram. 
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Figure 4: Variation of the critical temperature with q for some selected film 
thickness values of L = 3,4,6,10. The horizontal thinner line represents the 
critical temperature of the isotropic Heisenberg model for a simple cubic lattice. 



When the surface anisotropy is higher than that of the interior of the film, 
surface magnetization has a lower value than that of the interior. This situa- 
tion can be seen in Fig. ([Sj for a film that consist of XXZ-type surfaces on an 
isotropic interior, with L = 7. In Fig. ([5]), m stands for the total magnetiza- 
tion of the film, while mi denotes the magnetization of the surface layer and 
m4 denotes the magnetization of the middle layer. As seen in Fig. ([5]), rising 
q (which means rising anisotropy in the exchange interaction of the surfaces) 
depresses the surface layer magnetization to zero. But as seen in Fig. ([5]) (d), 
while the surface has zero magnetization at lower temperatures, rising temper- 
ature gives rise to non-zero magnetization value for the surface layer. Thermal 
agitations can give rise to surface magnetization which is different from zero at 
high anisotropy values in the exchange interaction of the surfaces. 



4.2.2 XXZ-type surfaces on a XXZ-type interior 

In this case, both of the exchange interactions of the system are anisotropic as 
XXZ-type. Here Eq. ([T5)) is valid but in contrast to previous subsection, r2 is 
different from unity. The problem of the behavior of the order parameter has 
four parameters, namely ri,r2,L and ksT/J. Thus, the critical temperature 
of the system depends on the number of layers of the system (L), anisotropy in 
the exchange interaction of the surface (ri) and the anisotropy in the exchange 
interaction of the inner layers (r2). 

The phase diagrams in a {kBTc/J, r2) plane can be seen in Fig. (J6l) for some 
selected values of g = r\lri- The thinner lines in each figure corresponds to the 
variation of the critical temperature of the corresponding bulk system which 
has spin-spin interaction (r-2)- The symmetry of this interaction is given in Eq. 
(fT5|) . Fig. ©(a) corresponds to the system that has Ising type surfaces on a 
XXZ-type interior, while Fig. ©(c) corresponds to the system with XXZ-type 
surfaces on a XXZ-type interior with the same magnitudes of anisotropy. The 
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Figure 5: Variation of the magnetization of the surface layer (nii), inner layer 
(7714) and total (m) magnetization with the temperature for L = 7 and selected 
values of (a) q = 2.0, (b) q = 3.4, (c) q = 4.5 and (d) q = 6.0. 



other two panels correspond to the system of XXZ-type surfaces on a XXZ-type 
interior, but in Fig. ©(b), the surface anisotropy value is lower than that of 
the interior {q = 0.5), whereas in Fig. ©(d), the surface anisotropy value is 
greater than that of the interior (g = 3.0). 

First, rising anisotropy in the exchange interaction of the interior (r2) de- 
creases the critical temperature of the system for any q, as expected. The 
important point is that the special point appears at the lower values of q and 
it disappears while q rises. Let us denote the coordinates of this special point 
as {ksT* /J,r2). The evolution of the special point coordinates with q can be 
seen in Fig. ^. As seen in Fig. ^, special point disappears just after the 
value of q = 0.83. Thus, if the ratio of the surface anisotropy to the interior 
anisotropy value is greater than q = 0.83, then thinner films have lower critical 
temperature for all values of the q. For the values that provide q < 0.83 special 
point can appear (as seen in Fig. (jH]) (a) and (b)) and in this case, for the 
anisotropy values that provide 7'2 > 7'2 with q < 0.83, it is possible to observe 
thinner films with higher critical temperature in comparison with that of the 
thicker ones. This special point moves on the lines can be seen in Fig. (O as g 
rises. 

We can see from Figs. ([6]) (c) and (d) that, if the special point is not present 
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in the system then ah the phase diagrams he below the phase diagram of the 
bulk system (denoted by thinner lines). Besides, after some value of anisotropy 
(g), the phase diagrams of the thinner films depress to zero, i.e. all transitions 
are of the second order (e.g. the curves labeled by A and B in © (d)). This 
is interesting since at a certain value of q, when the thickness of the film rises, 
second order transitions turn into first order transitions (compare the curves 
labeled by B and C in dH) (d)). 

We again note that, some of the phase diagrams in Fig. (0) depress to zero 
whereas some of them do not. Since the formulation can give only second order 
transition temperatures, the phase diagrams that does not depress to zero have 
the portions of first order transitions which we can not calculate within this 
formulation. 
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Figure 6: The phase diagrams of the thin film consisting of XXZ-type surfaces 
on a XXZ-type interior, in a (fc^Tc/J, r2) plane for some selected values of q 
and for the thickness values ofL = 3,4,6,10 



5 Conclusion 

In this work, the effect of the anisotropy in the exchange interaction on the 
critical properties of the thin films have been investigated. As a formulation, 
the differential operator technique and DA within the EFT-2 formulation has 
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Figure 7: Variation of the special point coordinates (a) rj and (b) ksT* / J with 
q for the thin film consisting of XXZ-type surfaces on a XXZ-type interior. 

been used. 

The system has two different exchange interactions and the problem has been 
handled by grouping the systems according to whether the surface/interior has 
anisotropic exchange interaction or not as: isotropic type surfaces on an isotropic 
type interior, XXZ-type surfaces on an isotropic type interior and XXZ-type 
surfaces on a XXZ-type interior. 

For the critical temperatures of the films, isotropic type surfaces on an 
isotropic type interior system exhibit special point which has coordinates (ksT* /J, rl)= 
(4.8910,1.325). This special point is related to the relation between the ex- 
change interaction of the surface layer and inner layers. Thus, the anisotropy 
ratio between the surface and inner layers have to effect on this coordinate. This 
effect has been shown in Fig. ([7]) for the system with XXZ-type surfaces on a 
XXZ-type interior, and we have concluded that after a ratio of q = 0.83, special 
point does not appear in the system. Besides, in the system with XXZ-type sur- 
faces on the isotropic interior, it has been shown that rising surface anisotropy 
decreases the critical temperature of the film and most affected films are the 
thinner ones. 

For the thcrmodynaniical properties of the film, in the isotropic case (i.e. 
isotropic surfaces on an isotropic interior), based on the ratio of the exchange 
interactions between the surface and interior of the film, surface magnetization 
can lie below or above the magnetization of the inner layers. This situation is 
related to whether that ratio is greater than the special point or not. For the 
values that provide ri > r* = 1.325, surface magnetization is greater than the 
inner layers and vice versa. For the XXZ-type surface with isotropic interior, 
surface magnetization is lower than the magnetization of the interior of the 
film for the values of q > 1.0, where q relates the anisotropy in the exchange 
interaction of the surface with the inner layers. It has been shown that when 
q increases, the magnetization of the surface layer depresses to zero at low 
temperatures, but due to rising thermal fiuctuation which comes from the rising 
temperature can create a non-zero magnetization. 
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We hope that the results obtained in this work may be beneficial form both 
theoretical and experimental point of view. 
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